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Abstract 



We study theoretically electron transients in semiconductor alloys excited 
by light pulses shorter than 100 femtoseconds and tuned above the absorption 
edge during and shortly after the pulse, when disorder scattering is dominant. 
We use non-equilibrium Green functions employing the field-dependent self- 
consistent Born approximation. The propagators and the particle correlation 
function are obtained by a direct numerical solution of the Dyson equations 
in differential form. For the purely elastic scattering in our model system 
the solution procedures for the retarded propagator and for the correlation 
function can be decoupled. The propagator is used as an input in calculating 
the correlation function. 

Numerical results combined with a cumulant expansion permit to separate 
in a consistent fashion the dark and the induced parts of the self-energy. The 
dark behavior reduces to propagation of strongly damped quasi-particles; the 
field induced self-energy leads to an additional time non-local coherence. 

The particle correlation function is formed by a coherent transient and an 
incoherent back-scattered component. The particle number is conserved only 
if the field induced coherence is fully incorporated. The transient polarization 
and the energy balance are also obtained and interpreted. 

78.47,+p, 71.23.-k, 72.10.Bg 
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I. INTRODUCTION 



The response of electrons in a semiconductor to a strong short optical pulse is an impor- 
tant case of transient non-equilibrium processes in many body systems. It is discussed in 
detail in monographs by Haug and Jauho |IJ and by Bonitz 0. Step by step, the technique of 
non-equilibrium Green's functions (NGF) has been recognized as the most reliable and phys- 
ically most transparent method for handling strongly non-equilibrium processes. However, 
because of the prohibitive requirements on computational means, its practical use was until 
recently connected with various approximation schemes reducing the full Kadanoff-Baym 
equations (KBE)or their equivalents to simplified quantum transport equations by means of 
various decouplings ("ansatzes"). The generic ansatz appears already in the classical works 
by Kadanoff and Baym 0] and by Keldysh M, and it was followed by a number of its re- 
finements. This work was very successful in describing a wide variety of experimental data 
almost quantitatively. Still, uncertainties concerning the validity and more subtle properties 
of the additional approximations remain. 

Since early times, however, there existed an important line of research ||, || on the direct 
solution of the KBE for transients, both in electron systems and in nuclear matter, where it 
was mandatory to use the NGF and related quantities as two time functions. More recently, 
efforts have concentrated precisely on a direct solution of the full NGF incorporating both 
the electron-phonon and electron-electron interactions, as discussed already in Ref. |2| and 
expounded in more detail in the proceedings [|7| . The need to work with two time functions 
limits the physical approximations for the self-energies, as analyzed in the general references 
given above. 

We will consider the third principal mechanism of electron scattering, caused by a static 
random potential. In random semiconductor alloys the strongest scattering is often caused 
just by the disorder, which then dominates the early stages of the transient response to a 
pulsed disturbance. The case of disorder has some interesting special properties while pro- 
viding an excellent model situation to study general aspects of non-linear transients [§] || 
10|1 ['!]• This is well illustrated on our two-band model semiconductor with disorder only 



acting in the conduction band [12]. A dipole optical transition transfers a valence electron 
with a sharp wave vector value k to the conduction band. For one alloy configuration, the 
conduction states are random and have no wave vector at all. In the one-electron picture, 
this is a fully unitary evolution in which the transition is made to many of the random states 
with a full preservation of the quantum coherence. Alternatively, a configuration average 



is performed [|T3| accounting for the fact that the macroscopic response of the electrons is 
self-averaging ||14|| . The static random fluctuations of the alloy potential play the role of the 
quantum field fluctuations of the many body theory ("frozen phonons"). Now, the transi- 
tions preserve the fc-vector, but the final states are renormalized and decaying, independently 
of the excitation process or of the electron distribution. In analogy with the phonon case, 
this may be termed a dark polaron effect. The optical jump injects an electron into a bare 
Bloch conduction state, and its transition into the dressed state (formation process) plays 
an important role for the short time dynamics of the coherent excitation. At longer times, 
the incoherent back-scattering sets on to compensate for the decay of the coherent popula- 
tion. All these features coincide with those occurring in a genuine quantum field case, and it 
appears that for the disorder case the full NGF treatment is the only reliable approach. In 



2 



paper ||, we have developed the basic theory for NGF in the presence of external fields and 
alloy scattering and solved the equations of motion in the coherent potential approximation 
(CPA, for the equilibrium case and linear response see Ref. ||13|| ) in a closed form. The solu- 
tion was limited to an unrealistic case of a sharp onset of steady illumination, however. In 
||, this solution was used for a numerical test of one of the successful ansatzes, the so-called 
GKBA fl5| . In an adiabatic approximation was developed for slowly varying pulses, 
employing the mixed time-energy representation, and in we were able to present early 
results for retarded propagators in the presence of pulses of a general shape. 

This paper extends our previous work to a full computation of the response to an arbi- 
trary pulse, although only in the weak scattering limit of the CPA, the self-consistent Born 
approximation (SCBA). 

The alloy scattering on a static potential is specific in that it is strictly elastic and instan- 
taneous, with no internal time structure. This has an important consequence. The electron 
propagators do not depend on the electron distribution, being basically a configuration av- 
eraged one-electron evolution operator. The propagators can thus be studied, analytically 
and numerically as well, in an independent and comparatively easy step. In the second step 
yielding the full NGF, the propagators serve as an input and it is sufficient to solve a single 
linear equation for, say, the particle correlation function, instead of the usual set of two 
coupled non-linear Kadanoff-Baym equations. 

On top of that, the use of this asymmetric set of Dyson equations |16| permits a direct 
interpretation of some features of the non-equilibrium process. By concentrating first on 
the propagators, we can study in detail the formation stage and the transition to the long 
time ("quasi-particle") regime. The self-energy can be separated into its dark part and 
the induced component, which reflects the coherent coupling between the excitation and 
the scattering. Proper inclusion of this coherence effect appears to be essential to keep the 
theory conserving and consistent. This becomes apparent in the behavior of the particle 
correlation function, in the loss of coherence in the non-equilibrium electron distribution 
and in the validity of conservation laws. 

We start, in Sec. [H]> by describing chemical disorder in semiconductor alloys ( II A ) which 
serves to justify the simple two band model with a short-range disorder of atomic levels 
acting only in the conduction band, as introduced in |IIB| . The effect of the alloy disorder 
in the dark is discussed in Sec. [TTI|. First, the SCBA is introduced and the renormalized 
bands are obtained ( [III A|) . For an actual choice of the model parameters, the renormalized 
bands, their damping, and the weak field one-photon resonance conditions are presented 
( [III B|) . The next subsection [III C| deals with the dark propagators in the time domain. The 
formation process leading from a bare electron to a quasi-particle is analyzed. The SCBA is 
compared with the exact short time cumulant expansion of the GF and the formation time 
is estimated. 

In Sec. [TV], the action of the light pulse is included. Specific properties of NGF in 
disordered systems are introduced in [TV A| . General relations for the Langreth-Wilkins form 
of NGF are in [IV A 1| . Next two subsections IV A 2 , IV A3| introduce the GF for a single alloy 
configuration and its configuration average. A closed form of the perturbation expansion is 
obtained. In |1V B| , the explicit equations of motion for all components of the NGF are given 
( [IV B 1|) . In |IVB 2| we separate the dark polaron effect from the light induced self-energy. 
Its properties are analyzed in [IV B 31 where the short time expansion is extended also to the 
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excited system. The meaning of the dark quasi-particle approximation is clarified in [IV B 4 . 
The actual computational schemes for both Dyson equations are presented in Subsec. [IV C| . 

Sec. |V| is devoted to numerical results. The actual choice of the pulse is specified (Subsec. 
V A| ). The retarded propagator is presented in [VB| . An example of the induced self-energy 



is shown in Subsec. |V B 1| and its role for the propagator is demonstrated in |V B 2| . The 
related coherence and time non-locality leads to a breakdown of the semi-group property of 
the propagator, Subsec. |V B 3| . The computed behavior of the particle correlation function 



is described in [VC| , presenting E < (t, t') in |VC 1| and an example of G < (k; t, t') itself in [VC 2 



Physically, we are particularly interested in the averaged one-particle density matrix, and 
in the observable averages computed with its aid. Sec. [VT] is divided into two subsections. 
In the first one, [VI A| , we consider the connection between the one electron density matrix 



and the particle correlation function, construction of the precursor kinetic equation and a 
proof of the particle number conservation in SCBA ([VI A 1| ) . Several observables and their 
average values are treated in VI A'2 . In the last subsection |V1B we show the results of 



the corresponding computations. First, individual matrix elements of g(t) are shown and 
discussed in |V1 B 1| . The last part |V1B 2| concerns observables. The light-disorder coherent 



coupling conditions the particle number conservation which hinges upon a full inclusion of 
the induced self-energy in the computation. The energy is not dissipated in the model, 
only transferred between the pulse and the electrons, while the dephasing of the interband 
offdiagonal components of the polarization need not be introduced by hand, as it appears 
to be linked with pronounced cancellations of the oscillatory integrand in the trace yielding 
the total polarization. 

Two appendices outline the steps necessary to incorporate into the theory random initial 
conditions (|A|) and observables represented by random operators (|B|). 

II. MODEL DESCRIPTION OF THE SYSTEM 



In this paper, we use the model semiconductor alloy introduced previously [§], ||10| . The 
simplifications inherent to the model bring it to a close parallel to two-band models with 
electron-phonon interaction Static alloy disorder in this picture may be understood in 
terms of "frozen phonons". The model belongs to the family of "independent band models" 
used to study optical transport in disordered systems [12| . 



A. Electrons in semiconductor alloys 

One aim of this paper is to turn attention to the importance of disorder effects in the 
optical response to short pulses; once the excitation will take place at regions of the band 
structure even moderately beyond the absorption edge, the alloy scattering in a mixed 
crystal will frequently become stronger than either the e — e or the e — p collisions and it 
will represent the principal mechanism of relaxation of the photoexcited population during 
the initial period. As an example, consider the zinc blende semiconductors. Even for a 
conservative choice, the Gai_ c Al c As system on the GaAs side, c < 0.1, pronounced alloy 
scattering effects are predicted |L7] . Stronger effects still can be expected for other systems, 



in particular some of the II- VI alloys . In all such materials, the basic alloying mechanism 
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is an isoelectronic substitution, by atoms from the same group of the periodic system. As 
a consequence, the random one-electron potential in the alloy is dominated by short range, 
basically site-localized fluctuations. In an orbital picture, this is best captured by a one- 
electron Hamiltonian in a localized orbital ("tight binding") representation which is 



by now standard, and considered quantitative |HJ, ||20fl . The tight binding Hamiltonian in 
the minimum basis of quasiatomic orbitals will typically involve (i) non-random off-diagonal 
hopping matrix elements responsible for bonding between the sp 3 hybrids, and for creation 
of the corresponding hybridization gap between the valence and the conduction bands, and 
(ii) random diagonal elements, quasiatomic levels q, with as many admissible values, as are 
the alloy components (see, for example, fll8"|)Jn the case of Gai_ c Al c As, the doping proceeds 



only on the cationic sublattice. In the Harrison parametrization [fLG 



e£?= - 11.37 eV, e£= - 10.11 eV 
e ga = _ 4.90 eV, eg= - 4.86 eV. 

If the alloy is random, that is, without any short or long range order, the random potential 
is fully characterized by the e values and the concentration c. In particular, the ratios 
(level difference) /bandwidth determine the alloy regime. If they are much less than 
one (as is the case for Gai_ c ALAs), a perturbation approach is justified. The ^-contribution 
to the level broadening is then proportional to the projected density of states of a given 
symmetry. In general, spectral weight of each orbital is divided between the conduction 
bands and the valence bands, so that, say, the disorder of the cation s levels affects both. 
The projected spectral weights are distributed unevenly, however. Around the gap, where 
the optical excitation is assumed to act, the conduction band is predominantly composed of 
cation states, and these mostly with the s-symmetry. The top of the valence band, on the 
contrary, is nearly anion p by nature. The cationic disorder, as a result, acts much more 
strongly at the bottom of the conduction band than at the top of the valence band. 



B. Two band model 

On the basis of these remarks, we will now develop a simplified model of the electronic 
structure of this system as used in Refs. ||, [ |TU| . 

1. We consider a two-band semiconductor with the gap between two isotropic parabolic 
band edges at the center of the Brillouin zone ("standard band structure") 

2. All many-body interactions are ignored 

In general, there are three mechanisms by which the disorder couples both bands: the c — v 
chemical hybridization mixing the conduction and the valence states, a statistical correlation 
between the electrons moving in both bands caused by the atomic disorder configuration 
they have in common, and a dynamical transfer of disorder between the bands due to the 
mixing of optically coupled valence and conduction states. We will eliminate both static 
coupling mechanisms by hypothesis. 

3. The disorder acts in each band separately, not mixing states of both bands ( "indepen- 
dent bands") 
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4. The disorder in the valence band is neglected ("cationic doping"). 

5. The effect of the optical pulse is restricted to a non-random interband dipole coupling 
treated in the rotating wave approximation (RWA). 

This model appears as reasonable for an optical excitation slightly beyond the absorption 
edge, but the band structure has to be extended to the whole Brillouin zone for two reasons: 
o the single site nature of the disorder spreads its effect homogeneously over the whole 
Brillouin zone, o a short strong optical pulse acts in an extended part of the Brillouin zone, 
(E—t uncertainty). This becomes apparent, for example, in the band renormalization caused 
by the disorder, and, in particular, in the short time quasi-particle formation process. The 
band parts remote from the Brillouin zone center can be described rather schematically, 
however, as only their gross properties enter the problem. 

Disorder does not mix the two bands. The band projectors of these "independent bands" 
are non-random and diagonal in both the Bloch and the Wannier (site) basis: 

P h = J2\bk)(bk\= £ \bi)(bi\ , b = c,v. (1) 

BZ lattice 

The full one-electron Hamiltonian for one configuration of an Ax_ c B c alloy has the following 
structure: 

H = W v + W c + V c + U(t). (2) 

The non-random configuration independent or averaged quantities are denoted by italics, 
while the configuration dependent ones by rounded (" \cal ") characters. 

Here, W v + W c is the Bloch Hamiltonian of the A crystal. Both band Hamiltonians are 
diagonal in the Bloch basis: 

W b = P b W b P b = J2\bk)e b (k)(bk\, b = c,v (3) 

k 

The random potential acts only in the conduction band and it is site diagonal: 

V c = P C V C P C = EM)e.H (4) 

i 

Here, are quasiatomic level shifts e A or e B . 

Finally, U(t) is a non-random inter-band dipole coupling of the pulse to the electrons, 
whose electric field is a linearly polarized harmonic wave with a basic frequency Q and an 
envelope E M $(t) with $ > 0, max$ = 1 ,E M is the peak value of the field: 

E(t) =E M cos(ftt)$(t)e (5) 

Assuming validity of the Rotating wave approximation (RWA), and making at the same 
time the approximation of a /^-independent transition matrix element, we get for U(t): 

U(t) = U cv {t) + U vc (t) = -e(r cv + r w )E(t) 

^ -$(*) (E \ck)Qe- int (vk\ + E \vk)Q*e im {ck\\ (6) 

I k k J 

t cv = P c tP v , etc., Q = |eE M (ck~ 0|(e, r)|uk~0) 
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The coupling strength parameter Q can be made real positive, Q = Q* > by rephasing. 
RWA allows to make a time-dependent unitary transformation 



= P C + P v e~ int (7) 
of the Hamiltonian (the Galitskii transformation |2~lf) to eliminate the rapid Q oscillations 



of U(t). Thus, Hit) H = OHO^ - Oihd t O*: 

H(t) = W v + W c + V c + U(t) 

e„(k) =e v {k) + Ml (8) 

U(t) = -Q$(t)(E|ck><uk| + Ekk)(ck|) 
Ik k ) 

The valence band is displaced by hQ upwards on the transformation. The interaction with 
light appears now as a transient hybridization of the transformed bands giving rise to the 
band splitting sometimes called the " Galitskii gap" . 

From now on, we will work in the Galitskii picture dropping the tilde accents over the 
operators, but we will keep the tildes to indicate the shifted valence band energies. 

Denoting the usual configuration average by (■■■), we may introduce the mean-field 
Hamiltonian H MF and the configuration dependent random field T>: 

H(t)=H MF (t)+V(t) (9) 
H MF (t) = (H(t)) (10) 

W)) = o (ii) 

In our model, U(t) is non-random, while V is time independent. Forming H MF then means 
that the virtual crystal part of V c responsible for the rigid band shift is transferred to the 
renormalized conduction band W c : 

H MF = W V +W c + U(t) (12) 

with 

e MP (k) 

W c = E\ck)e c (k) + (e)(ck\ 

k (13) 
V = E|ci)(ei_-_(e>),(d| 

di 

where (e> is the average atomic energy (1 — c)e + ce B . The A crystal will be doped with B 
atoms, thus we choose 



e 



^-0, e B = 6, (e) = c-5 (14) 



There are thus two alloy disorder parameters, concentration c and the level fluctuation S. 
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III. DARK ALLOY 



The disorder scattering is present even without the optical excitation, and the c-band 
states are correspondingly distorted. The excitation transfers electrons into these disorder 
dressed states, and the behavior of the light induced transition states can best be described 
differentially with respect to the dark dressed states. This situation is similar to the dark 
polaron effect in the initial state [pUfl . The poles of the renormalized dark alloy bands define 
the weak field resonance transitions, but for strong pulses, the short time deviations from 
the pole (quasi-particle) behavior are essential. 

In this section, we will consider this dark alloy case. This is an equilibrium situation char- 
acterized by a time independent 7i. While in general the electron response is conveniently 
expressed in terms of the non-equilibrium Green function, in equilibrium it is sufficient to 
consider the retarded propagator G R of one electron moving in the dark conduction band. 

The propagator then depends just on the time difference t — t', and may be analyzed 
without invoking the NGF formalism in the spectral representation introduced by Fourier 
transformation: 

G R (k, t)= f^r G(k, E + iO) • e- iEt ' h (15) 
J 2nn 



A. Dark alloy bands in SCBA 

A single conduction band with diagonal disorder is easy for treatment in any single site 



approximation, as described in detail in [13|. Here, we only sketch the Self-consistent Born 
approximation results. The Green function in the energy representation coincides with the 
configuration averaged resolvent, and is given by 

G(k, z) = —— (16) 

z - e k - (e> - E{z) 

Here, z denotes the complex energy, = e c (k) and the SCBA self-energy does not depend 
on k. It is determined from 

E{z) = 7 • F(z) 

F(z) = N- l YG{k,z) (17) 

k 

7 = c(l - c)5 2 

The Nordheim parameter 7 measures the scattering strength of the random potential, F 
is the local Green function(LGF); its imaginary part determines the conduction band DOS 
per cell by g c {E) = =|=7r~ 1 Im.F(.E' ± iO). Combining Eqs. (|16"D and fll7|) , we express this local 
GF in the alloy in terms of F a , the pure A crystal LGF, with a shifted argument: 

F{z) = F (z - <e) - E(z)) (18) 

The SCBA equation flTT| ) reads explictly: 
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E(z) = 7 • F (z - <e> - £{z)) 



(19) 



This can be solved by iteration. The dispersion law and the pole approximation, however, 
may be obtained explicitly. Let 8 grow from zero to its actual value. The band energies 
turn adiabatically into the complex resonance energies, = e c (k) i— ► z k = z c (k), given by 
the poles of the Green function ( |I6|) (on the analytic continuation to the non-physical sheet 
of the z-Riemann surface]: 



G-\k,z k ) 

Zk 





e k + (e) + E(z k ) 



This holds for any fc-independent self-energy. In the SCBA, 
pression for the SCBA dispersion law, 



Zk 



ek + <e> + 7 • F Q (e k + iO) 



(20) 

becomes an explicit ex- 
(21) 



The result looks like a most naive use of the non-selfconsistent Born approximation. 
Next, we rewrite Eq. ([16|) into a residuum form, 



G(k,z) 



Zk 



{z - z k ) 



R{k,z) 



(22) 



where Z k is the renormalization constant and R(k, z) is regular in the neighborhood of the 
pole. From ([19]) and (|20|), we get also the renormalization constant explicitly: 



1 + 7- ^F a (z) 



(23) 



B. Parametrization of the dark alloy model 

A sufficient input for the dark alloy includes o alloy parameters, 5 and c, o the c-band 
density of states per cell, g Q , from which the local GF F Q can be obtained by the usual 
spectral representation 

F (z)= [ d V ^— 9o ( V ) (24) 
J z — f] 

Because we will be interested mostly in the spectral region of the lower band edge, a model 
g a may be considered satisfactory, if it has the correct energy position of this edge and 
the curvature yielding the required (density) effective mass. Otherwise, it should have a 
proper bandwidth and it must satisfy the sum rule / dr]g (r]) = 1. One easy way, adopted 
presently, to adjust the DOS to these requirements is to write it as an expansion in terms 
of the Tschebyshev polynomials of the second kind. An example of F Q resulting from such 
procedure is shown in Fig.|l| Details are given in an appendix to ||1Q|| . The parameters 
used were the band edge e c (k = T) = Eq = 1.5 eV, the effective mass m* = 0.4 m e , the 
half-bandwidth w c = 6 eV. 
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The renormalized dispersion law Zk = z c (k) as a function of £/., is given by fl2ip . The 
renormalization has two parts: a rigid shift (e) = c5 of the whole band proportional to the 
alloy concentration, in agreement with the known notion of the rigid band model |T3| , and a 



complex component 7 • F a (ek + iO), whose real part gives the negative "polaron shift", and 
whose imaginary part is related with the electron life time by r « h/2\lmzk\. The polaron 
shift varies with concentration in the "Nordheim" fashion, proportionally to c(l — c). It is 
not entirely rigid, which leads to an energy and concentration dependent renormalization of 
the effective mass. This effect is minor in our case, however. 

In Fig. we display by thick lines the renormalized quasi-particle energies as a 
function of the bare band energy e&, taking 5 = 0.84 eV, and varying c in the range -j- 
0.5. The energy dependent band broadening is depicted in Fig. |] by lining each of the 
renormalized dispersion laws Kezk by thin lines at a distance ilmzfc. The rigid band energies 
are shown by dashed lines, and the polaron shift is given by the vertical distance between 
the full and the dashed lines. 

The corresponding renormalization constant is obtained from fl23|) and is shown in Fig|| 
This quantity is complex, and its module appears to be greater than 1, while the phase shift 
is negative. Thus, the long time behavior of the quasi-particle state will have the appearance 
of the simple Weisskopf-Wigner decaying state lagging somewhat behind in time. 

Returning to Fig. we discuss the plots (-.-.-) of the shifted valence band. With the 
basic frequency of the pulse specified, the valence band is shifted to e v (k) = e„(k) + HQ. Its 
crossing with the renormalized conduction band is the point of the weak field one-photon 
resonance (vertical transition): 

z x = z c (k x ) =e v (k x ) + fin (25) 

In the parabolic region, the two bare bands are related by 

C(k) = ^(^- ee (k)) (26) 

7Tl v 

(For simplicity, we will extend this as a definition for the entire model valence band.) Choos- 
ing m v = 0.6m e and taking first HQ = 2eV, we plot in Fig. the corresponding line. The 
crossing points with the renormalized c-bands define the electron and the hole energy at 
resonance, and also the respective excess energies easy to read as the separations between 
the crossing point and the band edges at the vertical axis. The bare energy of the crossing 
point serves to specify the position k x of the resonance in the Brillouin zone. 

The crossing points and the excess energies vary with the alloy concentration for a fixed 
Q, and it would be preferable to adjust the basic frequency used, so that either the crossing 
point k x or the conduction band excess energy be kept fixed. The vertical line laid through 
the c = crossing point defines, how the frequency should be varied with concentration in 
order to excite resonantly always the same part of the Brillouin zone. In practice, a nearly 
equivalent way of tuning the frequency Q — > Q x (c) may be to keep the excess energy equal 
to the difference hQ x (c) —Eg{c) constant, at least in the virtual crystal sense. The numbers 
corresponding to Fig. |2| are listed in Tab. | and give an estimate of the importance of the 
alloy effects: It is seen that even the rather moderate alloy parameters produce a strong 

effect, in particular for the electron lifetimes. For stronger alloy scattering, the quasi-particle 
concept becomes useless, because it would be damped before completion of its formation. 
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C. Time dependence of the dark Green function 



Inserting the pole representation (22) into the Eq. dl5|) , we have 

G R (k, t) = ^-Z k e~^ Zkt + ( ^- R(k, E + iO). e- iEt/h (27) 

III J ZTlh 

The first term prevails under typical conditions at sufficiently long times, and this corre- 
sponds to the quasi-particle behavior. 

To study arbitrary times, it is convenient to use an 'interaction picture' for G R by 
factorizing off the mean field (virtual crystal) propagator. For every k, we define (cf. (pf)) 

e MF = e k + <e> (28) 
G R (k, t) = j^e - * 6 - 1 x e"^ (i) (29) 

The additional factor e~I ? ^ describes the transformation to a quasi-particle state and its 
further evolution. We express it as a complex phase factor by introducing an "action" 
q = Re? + i Imq. Its real part describes the additional phase, the imaginary part measures 
the decay of the amplitude \G R \. This identification is verified in the long time limit. By 
(f27j), the leading terms in q represent the "polaron shift" z' k — e MF and the damping \z"\/%, 

q(t) « (4 - e MF )t + iz'i t + o(t) (30) 

For short times, we may employ the cumulant expansion of G, described in |[23|| , p. 3, ||24|| . 
Ch.2. that is to expand in powers of time: 

o- i<m = t^M p (31) 

p=0 v- 11 

^)=^Y, { —L^C P (32) 
p=i P- n 

Here, M p are moments of the spectral density of the state |k) calculated with respect to 
e MF , while C p are the corresponding cumulants. The moments, at least for small p, can be 
calculated directly without difficulty, proceeding, for example, similarly as in [25]. There 
are two options. Either the moments are calculated exactly, employing the fact that G R 
coincides with the configuration averaged evolution operator up to a 1/ih factor: 

G«(k,t) = i(k|(e4^>|k) (33) 
= T^e-l £MF * x (k|(e-H (w - eiIFM< )|k) (34) 

= -e-R*"' x K —^t mn - e MF l op ) p )|k) (35) 

w o p\n ' . ' 

M p 

Or the moments are determined within the SCBA. In that case, they should be identified 
with the coefficients of the Laurent expansion of G(z) at infinity, 
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G(k,z)=£-^- (36) 

p =0 \ Z fc MFj 



Now, we have to combine Eqs. (pf), ( |19|) and fl24T), using Laurent expansions of all quantities 
involved and determining the coefficients successively. 

Once the moments are determined in either way, the cumulants are obtained by com- 
paring both series in fl3"2"|). The results for p — -j- 4 are in the table [TT]; everything 
is expressed in terms of the moments of the VC density of states defined by fip°(E) = 
J dE g {E — (€))(E — E) p and of the alloy characteristics c a — cb, 8, 7- It should be noted 
that because of the choice E = e MF as the reference energy, the first moment and cumulant 
both vanish, Mi = C\ = 0. The vanishing first moment simplifies the cumulant expressions 
in all orders: Up to p = 2, the SCBA results and the exact ones are the same. Differences 
occur for higher p. The SCBA expressions use 7 as the only characteristic of disorder be- 
yond the mean field, as implied by (jlTI). In both approaches, the moments coincide with 
the cumulants up to the third order, and the first deviations occur in the fourth order as 
marked in boldface. 

We plot an example of the time evolution of the "phase" in Fig.[|. The gradual onset of 
^ for the shortest times justifies the conjecture that the one electron excitation starts as a 
bare (mean field) particle and develops the cloud of "frozen phonons" only in the course of 
time: 

G R (k,t) = I e -^(--^ c -* 2 +-) x e -W^ 2+ - (37) 

This behavior prevails only for very short times, and it is succeeded by a gradual transfor- 
mation to the quasi-particle mode, as seen from the figure. The characteristic time for this 
quasi-particle formation has the order Ti/{tk — Eg)- The energy entering this estimate is 
simply the energy distance from the nearest critical point of the density of states. In the 
present case, it yields a few, like 3, femtoseconds. Over this basic formative process, there is 
superimposed a much faster and rapidly damped evolution, with a characteristic time on the 
order of 10 _1 fs. We enclose these details into small windows which are then blown up as in- 
serts. The exact curves in the inserts are fitted by the lowest order cumulant approximations 
rather well. This testifies explicitly that the gross features of the whole band, as reflected in 
the lowest moments of the DOS are decisive here. In particular, it is interesting to inspect 
Fig. |]b depicting Im<^ related to \G R \. Because \Z\ > 1, the quasi-particle GF starts above 
the true one. Nevertheless, the latter quantity starts by a rapid decrease due to a spreading 
of the probability amplitude into the whole band. Only in the second stage of the formation, 
a recuperation takes place, and the exact Im^ starts to oscillate around the quasi-particle 
straight line. It is thus not easy in the present case to construct an interpolation form for 
the GF similar to that introduced in ||26|| . Instead, we proceed as follows. 



The averaged dark GF obeys the Dyson equation in the differential form, 

ih d t G R (t) = e MF G R (t) + J di£ R {t - t)G R {i) (38) 

Instead of this integral term involving the self-energy, we may define, in analogy to the 
Weisskopf-Wigner approximation, a multiplicative, single-time but A;— dependent, quantity 
0k(t) by 



12 



ih d t G R (t) = e MF G (t) + a k (t)G H (t) (39) 

An explicit "phase" integral equation relating the phase q with S R may be developed, but 
we will not pursue this task. 

Comparing with ([29]), we get the relation 

°k(t) = j f ,(t) (40) 



The last three equations combined will be used in Sec. |IV B 4] to separate the dark and the 
induced parts of the self-energy in an efficient manner. This will be based on the short 
time behavior of a^(t). As plotted in the other panels inserted in Figs. |]a, b, this effective 
self-energy varies strongly for very short times, but it tends quite rapidly, although in an 
oscillatory manner, to the constant value of the complex quasi-particle energy. 



IV. ILLUMINATED ALLOY 

The electron response to an incident light pulse is conveniently expressed in terms of the 
non-equilibrium Green function. The full formalism of the NGF is given in references [|lj or 
||. Here, only the directly important relations will be given. Our approach is characterized 
by two features. 

(i) We employ the NGF in the real time domain in the Langreth-Wilkins representation 
|i~6|l , that is we work with G R , G < , G A . We prefer the LW triplet over the usual Kadanoff- 
Baym pair G < , G > , because in the present case of a purely elastic scattering the dynamics of 
the propagators G R , G A becomes decoupled from the particle correlation function G < , and 
is much easier for numerical solution. The LW equations also have a slightly more "physical" 
appearance, suitable for interpretation. 

(ii) For non-interacting electrons, time evolution is strictly unitary even in a disordered 



system, while the genuine field theoretic description relates to configuration averages [JT3 
J27j. The non-equilibrium Green function needed here can be expressed in terms of config- 
uration averages in a direct fashion following ||. 

Below, we describe the LW formalism for the disordered alloys and the computational 
procedures for the electron propagator reacting to the light pulse, and for the particle cor- 
relation function. 



A. Non-equilibrium Green functions 



1. General relations and LW conventions 



We work in the real time domain and the NGF is given by a 2x2 matrix which 
we choose to be in the Langreth-Wilkins representation | l6fl : 



G 



G R G< 
G A 



(41) 
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Here, the component Green functions are one particle operators depending on two time 
arguments and satisfying the general relations 

G R (t,t') = {G A y(t',t) (42) 

G < (t,t') = -{G<y(t',t) (43) 

G>(t,t') = -{G>Y(t',t) (44) 

G R — G A = G > — G< (45) 

The symmetry relation ( f42"|) reduces the number of independent Green functions to two, 
while the last relation, Eq. ( pEo] ) {spectral identity), links the LW representation with the 
usual Kadanoff-Baym pair G < , G > . 

The so-called LW algebra for the individual components of (|4ll) follows easily from the 
usual rules for matrix operations, For example, introducing also the self-energy matrix S, 
we have 



Thus, 



EG 



{£ G} R 
{SG}< 







G R G< 


S A 




G A 



E R G R 

E R G < 



E<G A 



(46) 



(47) 
(48) 



We use the convention that all quantities with no time variables shown are double time. 
Thus, S E a ,p = £ X (t,t'), H UF # H MF (t) 5(t — t f ), etc. Multiplications mean an operator 
multiplication and a time integral, like E R G R ^ / diS R (t,t) G R (i,t'). 



2. Single alloy configuration 



For non-interacting electrons, time evolution is strictly unitary even in a disordered 
system: Each alloy configuration gives rise to a random Hermitian one-electron Hamiltonian 
Tl(t), which depends on time due to the effect of external fields (light pulse). The non- 
random configuration independent or averaged quantities are denoted by italics, while the 
configuration dependent quantities are denoted by rounded characters. 

For one configuration, we introduce the usual evolution operator S(t,t') corresponding 
to H(t), and the one electron density matrix g(t). Selecting the earliest admissible time to, 
we introduce p(to) as the initial condition for the NGF. Then, for all times t, t' > to, the 
random single configuration NGF is 



g(t,t') 



^ S(t, t')$(t -f) - i- S{t, t ) g{t ) S(t , t') 
-^ S {t,t')${t>-t) 



(49) 



Two properties of Q are very important: First, the propagators Q R , Q A do not depend 
on the particle distribution, and, for propagation from t' to t, they only depend on 7Y(t) 
with t between the two times. Second, the correlation function Q < does in fact not depend 
on the choice of the "initial time" t , if the relation g(t) = iS(t,t ) £>(to) <S(to,t) is employed. 
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3. Configuration average of Q 



Observable quantities are presumably given by the configuration average of the particle 
distribution function p = (g)( see Sec. |VI], however); to obtain this quantity, we introduce 
the averaged non- equilibrium Green function starting from the Eq. ([49|) : 



G = (Q) 



(50) 



These configuration averages ought to be treated by field theoretical means [13|], |p? 
The two propagator components are simple, for X = R, A, 



1 



G x (t,t') = (g x (t,t')} = ±-(S(t 1 t'))$(±[t-t'}) 

Now, we obtain a semi-explicit expression for 

g< = (g<) = ihig*^) (t > t , t' > t ) 



(51) 



(52) 



Here enters the initial condition at t = to; the initial one-electron distribution is represented 
in the double time form as 



oPo = p(t )5(t - t )S(t' - t ) (53) 

This initial condition can be random, that is configuration dependent. We will exclude this 
possibility in our model by hypothesis. In this paper, we consider the case that the initial 
distribution p(to) is non-random. By this, we eliminate the problem of initial correlations 
p8| . For the non-random initial condition, the averaged product in (|5^) can be written as 

G< = ihG R oPo G A + G R E<G A (54) 

COHERENT INCOHERENT 

This equation defines E < , which is a regular two-time kernel in this case. The Eq. (|54"|) 
is then nothing else than the <-component of the Dyson equation for the non-equilibrium 
Green function matrix fl41|). It is a good example of the suggestive LW form of the theory. 
It has the full R <-> A symmetry and a causal structure. The initial state evolves coherently 
in resemblance to unaveraged g < (cf. Eq. (|52"D), but it is evanescent. The population is 
replenished by an incoherent backscattering due to S < . Eq. ( |54"D is also a starting point for 
developing equations of the Bethe-Salpeter type and the quantum transport equations. 

Given these results, we introduce the corresponding Dyson equation in the differential 
form by starting from the equations for the unaveraged Green function and performing the 
configuration average term by term. 

The random Hamiltonian can be divided into its configuration average, the mean field 
Hamiltonian, and the fluctuation potential, that is the configuration dependent random field: 
'Hit) = H MF (t)+U(t), as introduced by Eq. (|9|). Both parts of Hit) might be time dependent 
at this level of generality, in contrast to the specific model, Eqs. (0), fll3| ) which we employ 
for explicit calculations in this paper. 

By the Eqs. (0), (0), and (|), 
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(ihdt - H MF )g - x>g = i (55) 

(ihdt - H MF )G - £G = I (56) 

Here, d t # ^5(t — t'), I# t^I # 5 aij3 l OP 5(t - f), and £> # X>(t)J. On going from (||) 
to (|56|), we define the self-energy matrix by setting (T>Q) e EG. The possibility of such 
definition is itself another form of saying that the initial condition is non-random. A brief 
discussion of the random correlated initial condition is deferred to the Appendix |A|. 

The Dyson equation QSBp has to be supplemented by an actual form of the self-energy in 
order to become closed. In this paper, we will use the Self-consistent Born approximation 
(SCBA), which is a generic self-consistent conserving approximation appearing as a lead- 
ing term either in the self-consistent multiple scattering expansion, or in the self-consistent 
perturbation expansion. The relationship between SCBA and the Coherent potential ap- 
proximation within the multiple scattering context is discussed in M. Here, we use the 
self-consistent perturbation approach. This is based on the integral equations 

£ = G + G (£>-£)£ (57) 
£ = G + £(£>-£)G (58) 

which follow from Eqs. fl55|) , (|56|) , and which can be iterated to yield a renormalized Born 
series in terms of the full propagators Q and the renormalized perturbation T> — S. To 
obtain a closed equation for the self-energy, we iterate (|57|) once using (|58D and configuration 
average respecting Eq. ([IT]): 

£=(£>££>) - £G£ (59) 
This yields an iterative expansion in terms of T> and Q 



S = (X>GX>)+(T>GX>GX>)+ ( 60 ) 
+ (T>GT>GT>GT>) - (X>GX>)G(X>GX>) - (T>G(DGT>)GV) + ... 

The expansion starts from the second order in V, which is our SCBA. In the fourth order, the 
expansion already deviates from a simple geometrical series. In the language of diagrams, 
the first subtraction excludes a reducible diagram, the second suppresses double counting of 



self-energy insertions in the inner propagator lines [|l]|. The physical approximation for X, 
like a termination of the expansion at some order, is made consistently for all components 
of the self-energy, which guarantees its self-consistent conserving property. This is also the 
case of the SCBA. 



B. Working form of the Dyson equations 

Now we write explicitly the equations for all components of the NGF in the form suited 
for numerical solution. These are obtained as the corresponding components of the matrix 
equations Q5"6D and flBTjp with the use of the LW rules ( fl6"D - (fjE§|). For each component, we 
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have to specify: o the equation of motion; o initial conditions; o the specific form of the 
SCBA self-energy. 

We repeat that the equations will be written in the Galitskii picture with the convention 
that 

the operators, Green functions, etc., will be written without tildes. To return to the 
Schrodinger picture, we will use the subscript . . 

1. Explicit equations for G R , G < 

For the propagators, we get 

(ihd t - H MF ) G R > A - E R ' A G R ' A = I (61) 

G R (t = t' + 0, t') = ih~ 1 .l 0P (62) 

Eg A = (VG R,A V) (63) 

In fact, only one of the propagators has to be found directly, as the other one is given by 
©■ 

For the particle correlation function, the three relations are as follows 
The equation of motion has no 5- singularity: 

(ihdt - H MF )G K — E R G < + E<G A (64) 

o It has to be integrated for all times t > to for any fixed t' > to- Thus, the initial value is 
G < (t ,t'). For a non-random initial distribution p(to), we get immediately 

G<{t ,t') = p(t )G A (t ,t') (65) 

O The SCBA self-energy E < has a form analogous to the propagator component: 

r< BA = (VG<V) (66) 

We will return to the equations flBTD - ( |66"D below and show the methods of their practical 
handling and the computed electron dynamics. Before that, we will make our model fully 
explicit . 

First, we select the initial distribution function by assuming that the valence band is 
completely occupied, the conduction band entirely empty: 

p(t )= £ \vk){vk\ = P v (67) 
ke bz 

As indicated, the p(t ) operator coincides with the projector on the valence band states, 
which is non-random and alloy composition independent by assumption, so that (B3) applies. 

Second, all equations have to be specialized for the two band model of Sec. |II B| . Thus, 
all GF, etc., will be fc-diagonal 2x2 matrices: G x ^ \\G^ b (k;t,t')\\, {X = R,A,<; a,b = 
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c,v; k e BZ} and similarly for the self-energy. Disorder is acting only in the conduction 
band, so that only the cc components of the self-energy are non-zero. The self-energy 
becomes fc-independent in the single-site approximation. By this, the general relation fl6"5|) 
simplifies to a scalar equation in the spirit of the dark relations (1171) : 

F c x c (t,f) = iV- 1 E GX(k;t,t') (68) 
keBZ 

7 = c(l-c)<5 2 



2. Separation of the dark self-energy 

In principle, we have to solve the Dyson equations (|6lD for G R and fl64[) for G < , and to 
simultaneously obtain the self-energies from the self-consistent relations fl6"3|), (|66|). These 
equations will first be modified, however, by singling out the dark constituent of the self- 
energy. For these modified equations, a self-consistent cycle leads to their stable solution. 
The self-energy will be divided into its dark and induced components not only for com- 
putational, but also for physical reasons. The dark component analyzed in the preceding 
section describes the polaron effects including the short time dynamics and the long time 
quasi-particle behavior. The induced component reflects the coherent coupling between the 
excitation and the scattering. Proper inclusion of this coherence effect appears to be essential 
to keep the theory conserving and consistent, as is documented below. 

We will separate out the dark part of the retarded self-energy; what remains, is the 
induced part: 

yR _ yR i yR (rq\ 

^ DARK ' INDUCED l UlJ V 

No such separation is needed for the particle component, because all of it is induced in our 
case: 

= ^INDUCED (70) 

For convenience, the subscripts will be shortened to D and I. 

The equations (0) and ( |6^D for the retarded propagator become: 

Hqp 

ihd t G R - '(H MF + G R = Z R G R , t>f (71) 
E R = (V(G R - G R )V) 

The dark self-energy is known once for ever for a given system, and it appears as a permanent 
renormalisation ("polaron shift") at the l.h.s. of the Dyson equation. There is no coupling 
between this renormalizaton and the light. Without E R , the pulse would act as if it were 
exciting preexisting dark quasi-particles. Furthermore, the dark self-energy acts on the 
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propagator almost as if it were local in time. For all these reasons, we denote the operator 
on the lhs of the Dyson equation by H QP ... an effective quasi-particle Hamiltonian. 
Only the Dyson equation ( |64D is modified for G < in view of Eq. (|70|), 

Hqp 

■ * (72) 

ihd t G< - (H M¥ + E R ) G< = r« D G< + E<G A t > t Q 

while the equations (|66|) , (|65f) remain without change. This triplet Dyson equation-SCBA- 
initial condition is similar in many respects to the system for the retarded component. In 
particular, the lhs of the Dyson equation has again the "polaron" interpretation of dark 
quasiparticles excited by the light pulse. 



3. Induced part of the retarded self-energy 

To see in better detail the structure of the induced part of the self-energy, we rewrite it 
explicitly: 

Z?Jt, = TiV- 1 £ (c£(k; t, t') - G« D (k; t, f )) (73) 
keBZ 

Three points are apparent. 

First, E R (t,t f ) = 0, if the interval [£,£'] has no overlap with the time span of the pulse, 
because then G R = G R . Thus S R vanishes in the quiescent time domain, and need not be 
computed for these times. This SCBA result coincides with its exact counterpart. Namely, 
in an exact treatment, the Green function is essentially a configuration average of a unitary 
evolution operator. The evolution is governed by a time local Hamiltonian, and the con- 
figuration average is also instantaneous. As a result, the averaged propagation from if to t 
cannot depend on external fields acting outside this interval, just like in the approximate 
treatment. 

Second, the excitation by a pulse is not homogeneous in the whole Brillouin zone, acting 
predominantly in a strip of a width ~ Q around the resonant energy. In the parts of the BZ 
far from this strip, the excitation is weak, and the full Green function deviates but weakly 
from the dark one. This averts the necessity to integrate over the whole zone - in contrast 
to the dark self-energy itself. Thus, only the physically relevant k— vectors are involved, and 
the details of the distant parts of the band structure have only a minor importance. 

This, in turn, reduces further the time domain, where the induced self-energy deviates 
from zero. While it vanishes exactly for initial times t' after the pulse, we may conclude 
now that it is also negligible for t' well before the pulse. Namely, both G R and G R are 
exponentially damped at a rate roughly corresponding to the resonant energy for all relevant 
k; the whole induced self-energy is thus exponentially small in t — t' as t reaches the time of 
the pulse. 

Third, the induced part S R is continuous for equal times, t = t'. Each of the self-energies 
has aj/ih jump at equal times, and these two jumps exactly compensate. By extending the 
moment analysis of Sec. | [I G| also to the full G R , it is possible to improve the short time 
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estimate of E^. The moments depend on the inital time t', and we quote without details 
that Eq. ( j37)) is generalized to 

(i-iT + gW)) 1 ^^----) ( 74 ) 

1 i eMF(k)( ^ )x 

in 

For the shortest times, the contributions to the "action" due to light and to the disorder 
are independent and additive. Introducing this expansion into the Eq. (J75D, we deduce the 
short time behavior of E^- in the form 

£2A = 7Q W)K^ } M F„(t - O^f^ + ©((t - t') 3 ) (75) 

Here, F a (t) is the Fourier transform of the pure crystal local Green function, see flTB|). This is 
an important result, because it shows that the complicated short time behavior of the dark 
Green function does not enter the induced part of the self-energy, which gives it a welcome 
robustness with respect to approximations for Green functions. 



4- Dark quasi-particles 



So far, the definition H MF + El 



H QP had only a symbolic meaning, and was not 



associated with an explicit introduction of the quasi-particle picture into the computations. 
There are reasons, both practical and theoretical, to investigate this possibility in more 
detail. On the practical side, the dark self-energy varies strongly in a very short time 
interval. To include it properly means to perform the integrations on the lhs of the Dyson 
equation with a great care. At the same time, it is clear that beyond a short formation 
time, the dark propagators will assume the Weisskopf-Wigner form of a renormalized time 
exponential. Instead of directly using these WW propagators as a basis, we recall their phase 
form discussed in Sec. [Ill C| . The GF of that section coincides with G^ c , the time t is replaced 
by t — t' in the non-stationary situation in the Eq. (P§|). For large time differences, the 
action of E^ will approach the WW limit as given by (plD, — > Ofc(oo) = 7 • F Q (ek ± iO) 
regardless of the presence of the external field. On the other hand, for short times t — t', we 
have shown in Eq. ( [74]) that the phase variations caused by disorder and by the light are 
additive. Thus, the effect of the disorder and of the external disturbance appear as additive 
both in the short time and in the long time limits. Interpolating this behavior, we arrive at 
the following quasi-particle approximation: 



H, 



QP\k 



e c (k) + <e) + a k (t 

-QHt) 



t') -Q*(t) 
e v (k) + hn 



(76) 



It is understood that this effective Hamiltonian acts on a Green function with time arguments 
t,t' and crfc(t — t') acts as a multiplicative quantity. If the renormalization effects are not 
strong (slow energy variation of Fq), a further approximation, pure WW o~k(t — t') ps (00), 



may be tested, as we did in [O. 
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C. Computational schemes 



1. Computational scheme for the propagator 

We recall once more the important point that the propagator dynamics is fully inde- 
pendent of the particle distribution in our model with elastic scattering. The propaga- 
tors and the related self-energies can thus be computed from Eq. ([H]) in an independent 
round of the whole calculation. This is a comparatively less demanding computational 
task. Only the retarded GF and self-energy need to be obtained in a self-consistent cy- 
cle. The results can be stored, the advanced counterparts are obtained using the symmetry 
(0): G A (t,t') = [G R (t',t)]\ E A (t,t') = [E R (t',t)]l In Fig. §, we show the flowchart for 



computing E R . For solving the Dyson equation, we used a fourth order adaptive Runge- 
Kutta-Fehlenberg solver. Interesting is the seemingly counterintuitive order of time loops. 
This is dictated by the fact that we are bound to make one step in the floating time, 
t — > t + At, for all intermediate initial times t > i > t', so as to be able to calculate the 
product E R G R — ► / diE R (t + At, t) G R (i, t') for the new time t + At. 



2. Computational scheme for the particle function 



On the whole, the Eqs. and ( 72]) have an analogous structure. There are also 



marked differences. For G K , the integration of the Dyson equation starts at t = to for any 
fixed t' . The advanced propagator enters ( [72] ) at two places, in the initial condition and on 
the rhs of the Dyson equation in the integral / diE < {t, t) G A (i, if). Thus, G A (i, t') is needed 
for to < t < t'). It is readily available from the retarded propagator G R (t',i) calculated and 
stored beforehand, at the retarded stage, by means of the crossing symmetry (0). Another 
difference is that the self-consistent equations for G < and E < are linear, once the propagators 
are known, so that no iteration for obtaining E < is needed in principle, and the integration 
process runs only once. An iteration of the process appears as necessary on the practical 
level, however. Namely, in contrast to the retarded case, where E R (t,t) = is known, there 
is no universal value of E < (t,t), and this quantity must be estimated by extrapolation and 
refined by iteration. The whole process is outlined in a flow chart shown in Fig. []. For the 
particular Hamiltonian (§) and T> given by (|13D , the explicit form of the integral for E K is 



E<(t,t')= 1 N- 1 G<{W) (77) 

kGBZ 

7 = c(l - c)5 2 

The number of primitive cells is denoted by N, the Nordheim parameter 7 as the second 
cumulant of T> measures the disorder scattering strength in the SCBA. No subtraction of 
the dark GF is needed, because without illumination the electrons cannot leave the valence 
band, G< = P V G<P V , while disorder is confined to the conduction band. All of the cc- 
component of the full G < is induced, as stated in Eq. ( ffO|) . The integration in ([77]) extends 
in fact only over a narrow shell around the one-photon resonance, just like in the integral 
([HI) for E R Q . There, however, the subtraction is crucial. 
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V. NUMERICAL EXAMPLE 



In this section, we show a representative numerical example for the alloy model defined 
in Sections [II B| and |IIIB| . The bare A-crystal band structure is characterized by the band 



gap E, 



G 



1.5 eV, the electron and hole masses m r = 0.4 to* 



m, ; 



0.6 m P . The bands are 



"proportional", with band widths 12 eV, and 8eV, respectively. The band density of states 
is given in [[IIB} The alloy parameters are: c = 0.05, 5 = — 0.84eV, so that (e) = 0.042 eV 
and 7 = 0.034 eV 2 . We use the basic frequency TiVL = 2eV for a pure A-crystal and adjust 
the frequency for the actual c to 1.75 eV according to Table [I] The exciting light pulse is 
specified in |V A| below. 

All parameters are selected so that they lead to pronounced but moderate effects. The 
formation times are shorter than either the pulse duration or the electron relaxation time, 
and the photoexcited electron density is sufficiently low that the e — e relaxation time is 
longer than both these times. 



A. Pulse parameters 

We have to specify the pulse strength Q and its shape 3>(i), see (§). The pulse strength 
is Q = |eE M (ck ~ 0|(e, i)\vk ~ 0) by (H). The x-matrix element is estimated as the lattice 
spacing ao ~ 10 a.u., corresponding to the momentum matrix element of about 0.5 a.u. As 
max$ = 1, we get the estimate E M [Vm _1 ] = 1.89 x 10 9 Q[eV] for the peak value of the 
electric field. 

We use the "sech 2 " shape of the pulse and control its duration by a single parameter t P : 

$(t) = ¥#-ti/ 2 )/tp) 
<p(t) = 4/(e' + e^) 2 

The ti/2 shift is defined by the condition $(0) = |. This shift and the full width at the half 
maximum of intensity are related to t P by p9[ 



t 1/2 = t P ■ 0.88137 
^fwhm — ^"P ' 1.21169 
A convenient integral measure of the pulse is the resonant Rabi phase, 

+ 0O 

<p K = J dtn K (t) (so) 

— oo 

where the Rabi frequency is introduced by 

±M\(t) = Q<f>(t) (81) 



For our pulse shape ( [78] ) this gives 

<p R = 4hrH P Q 



(82) 

3.3012^. ^fwhmQ 



We use a pulse with characteristics according to the Table fTJ. 
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B. Propagators 



1. Calculated E^(t,t') 

G R and E R are given by the Eqs. (71 ) and (|73|). These were solved as explained in 



IV C . First, we present in Fig. |7| the imaginary (dominant) part of the induced self-energy 
as a function of both times t, if . This calculated shape manifests the properties predicted in 
Subsec. |1V B 3| . As a function of t for a fixed initial time t', it is strictly zero for t < t' and 
it continuously assumes non-negative values on crossing the diagonal t = t' . For t > t', it 
shows a marked memory effect, that is time non-locality, or temporal coherence, and it ends 
with a tail corresponding to the dark decay rate of the excitation after the end of the pulse. 
In the perpendicular direction, as a function of t' for t fixed, the ImZ^i profile is symmetric, 
decaying exponentially for times t' well before the pulse arrival and ending continuously but 
abruptly at the upper limit. Across the time diagonal, the damping is much faster, as seen 
both in the main plot, and in the rotated insert. To motivate this, we note that in the 
Wigner coordinate t — t', the exponential damping has a doubled decrement: G R, s in ( |73|) 
depend only on t — t' once both times are large enough, and for t + 1' — const, we have 
t — t' = It — const. Along the time diagonal, as a function of t + t', ImE R repeats well the 
shape of the pulse. 

The induced self-energy has a much smaller magnitude than its dark counterpart, but 
it persists much longer. As a consequence, it is difficult to qualitatively guess the relative 
importance of both components of the self-energy for the resulting propagators. 



2. Calculated G R (t, t') 

In the two-band model, the propagators are given by 2 x 2 complex matrices 
||G„ b (k;t, We show the behavior of the GF for just one point in the BZ, namely 
k m k x , by examining the time evolution of two quantities, \G R C \ 2 and |G^,| 2 . For a time 
development which started from a pure k x c- or from a pure k x t>- state, respectively, they 
measure the probability that the particle stays in its initial state. In Fig. ||, we plot the 
two quantities as a function of t and t'. For equal times, they start from a constant value 
h~ . For initial times well out of the pulse time interval, the decay of the c state resembles 
the dark alloy, while the v state does not decay at all. Within the pulse region, the very 
earliest evolution is still dominated by the dark self-energy. Soon, however, the external 
field controls the evolution, and the Rabi oscillation develops. 

This qualitative picture does not change, if the induced part of the self-energy is switched 
off, as depicted by plotting two surfaces in each panel of the figure, which resemble each 
other quite closely. The differences are but quantitative. In other words, the dark quasi- 
particle approximation of Sec. |IVB4| neglecting S R is qualitatively satisfactory. It has its 



serious problems, however, as we show now. 
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3. Testing the semi- group property 



There exists a sensitive test of the quasi-particle behavior of the retarded Green function, 
which leads to an easy and convincing proof of the importance of the induced self-energy 
||1 1|| . The true Weisskopf-Wigner dark quasi-particles, with zero illumination, would have 
an exponentially damped modulus. A generalization proper for the dark quasi-particle 



approximation in the WW limit, as defined in Sec. [I V B 4| , employs the time-locality of 
the WW quasi-particle Hamiltonian to justify the semi-group multiplicative property of the 
propagator. This is summarized in the following equations: 

ihg R (t,t r ) = ihg R {t,t")-ihg R (t",t') 
ihG R (t,t') ihG R (t,t")-ihG R (t",t') 
ihG R w {t, t') = ihG R w (t, t") ■ ihG R w {f, f ) (83) 

t > t" > t 1 



For the unaveraged GF, the equality follows from Eq. (f±9|), while for the exact averaged GF, 
the inequality expresses the fact that a product of averages does not equal to the average 
of a product of two quantities. For the quasi-particle case, the equality is restored by the 
time locality of the effective Hamiltonian (|76|). We study again G R for the resonant fc-state 
starting this time at t' = — 0.2ps. It is preferable to examine the squares of the full coherent 
amplitudes, \G R C \ 2 + \G R C \ 2 and |G\^| 2 + |G^| 2 . In Fig. |9|, we see how the c-band state 
initially decays with the dark decay rate, but later one Rabi flip is apparent. The f-band 
state is stationary at first. The pulse leads to its substantial depletion; a hint of the Rabi 
oscillation is also noticeable. The left panels incorporating U R and the right panels using 
the WW quasi-particle approximation lead to a picture qualitatively resembling the exact 
result to a large degree. 

Now we start testing the semi-group property (|83|). As marked by squares on the profile 
of the pulse, several instants t" were selected, at which we factorized the GF and plotted in 
various thin lines the time evolution corresponding to ihG R (t,t")ihG R (t" , —0.2) for t > t" . 
While these thin lines coalesce with the full line of the non-factorized Green function in 
the quasi-particle case, the factorization for the true GF appears as approximately correct 
only for factorization times well before or well after the maximum of the pulse, when, of 
course, we deal with a nearly dark evolution, so that the time non-locality of the self-energy 
does not play an important role. We will see in Sec. |V1B 2j , how these deviations from the 
quasi-particle behavior influence also the evolution of the photo-excited electron population. 

It should be pointed out that these results were obtained in the WW limit, because 
the discrepancy between both cases is particularly clear. Similar differences caused by the 
neglect of E R appear also in the general case and are not peculiar to the WW approximation. 



C. Particle correlation function 
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1. Calculated E<(t,t') 



The 'less' component of the self-energy is the central quantity for a full description of the 
excitation and transport processes generated by the pulse: in the Bethe-Salpeter equation 
E < appears as an irreducible vertex incorporating all correlations, which are in the present 
case due to the configuration average of the coherent multiple scattering of an e-h pair out of 
the equilibrium. As stressed already in Sees. [IV B 2| and [IV C 2| , E< would be zero without 
an excitation to the c band which activates the disorder scattering. E < has a coherent 
component resembling the pulse shape and width along the time diagonal t — t' = 0, as 
shown in Fig. |TD|. This is similar to r^ D , Fig. [7]. A marked new feature, however, is the 



appearance of a tail of E < which persists beyond the pulse duration. This is connected in a 
self-consistent manner with the photoexcited population in the c band. For the SCBA, we 
may see it clearly from the Eq. (|77|), by which E K is proportional to TrGf c . 

Finally, we look in the direction of the time diagonal (the arrow in the figure) to inspect 
the profile of S < along the other, t + t', diagonal. In the small inserts, the plots of the real 
and the imaginary part of the self-energy are correspondingly rotated and the symmetry 
relations equivalent to (f43|) are demonstrated. This represents a non-trivial check of the 
numerical approach used, because the integration procedure based on the differential Dyson 
equation ([72] ) is not a priori symmetric with respect to both times. 



2. Calculated G<(t,t') 

Fig. |y] shows just one representative Gf c (k; t, t'), namely that for k near the one-photon 
resonance. Symmetries ( f43|) are verified as a feature of individual terms of (|77J), as shown 
in the format of Fig. [H| We will concentrate on the even part hlmG^ c (k;t,t'), because its 
time diagonal equals to the induced population of the |ck) state. The induced transient is 
followed by pronounced structures slowly decaying along both time axes t, t' and by a slow 
rise tending to saturation along the time diagonal t = t' . A closer look shows about one full 
wave along the diagonal, in agreement with the Rabi phase of our pulse being adjusted to 
about 2tt. The side wings display only one rise followed by a decrease suggesting a phase 
variation ir. To understand the origin of this behavior better, we split in Fig. [12| the 
Green function into its coherent and incoherent parts in accordance with (|54]). The side 
arms of lmGf c appear to originate from its coherent part ImGp 0H . In fact, for one time 
fixed, the profile along the other time axis is proportional to G R , so that its oscillation is 
given by the Rabi phase halved, and the attenuation at long times is h/lm U^ ARK . Along 
the time diagonal, G< 0H is bilinear in G R , and the attenuation time equals now to the dark 
quasiparticle lifetime h/2lm S R ARK leading to a much faster decay. The incoherent part 
lmG< C0li appears with a time lag due to the multiple scattering delays. It is concentrated 
to a zone along the time diagonal, where it gradually supersedes the coherent part, so that 
it is responsible for the final rise of the full ImG < in Fig. [IT[ . 

This description is somewhat oversimplified and would be correct only in the dark quasi- 
particle approximation of [IV B 4| , in which the exponentially decaying quasiparticles are 



driven by the incident light. In fact, the Eq. fl72[) goes beyond the quasiparticle picture 



by incorporating coherence between the light and the disorder as captured by E R D . This 
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coherence was shown to be essential for a proper description of the propagators. It is equally 
important in the Eq. ( f72|) for G K . We compare for clarity the coherent parts ImG< 0H with 
and without including XJ^ D in Fig. O. While the two surfaces are superficially similar, 



they differ markedly at long times. It should be noted that while E^ D itself is a transient 



quantity, its coherent effects persist. This is well illustrated in Fig. [H] demonstrating the 
course of the full ImG < along the time diagonal. When ([72]) is solved without X^ D , the 
asymptotic value of the state occupancy differs from the full solution. This has important 
implications for the particle number conservation, as will be discussed in detail in the next 
section. 



VI. PHYSICAL PROPERTIES OF THE SYSTEM 



In general, observable physical properties are given by the equal-time limit of the particle 
correlation function In the disorder case, this correspondence is not so simple in case 
the observable itself is represented by a random operator. For our model, however, the most 
important observables can be reduced to traces involving non-random operators. 



A. One-electron distribution and observables 



For each configuration, the one-electron density matrix at a time t is given by g(t) = 
S(t,t )g(t )S(to,t). Thus, we have 

g (t) = -ihg < (t,t) (84) 

pit) = (g(t)) = -ihG < (t,t) (85) 

Consider now a one-electron observable represented by an operator X, which must be taken 
as configuration dependent in general. Its mean value susceptible to measurement is given 
by a double average, configuration and quantum statistical. It must in general be written 

as 



«*» = Tr(*£> (86) 

Thus, the configuration average concerns the product Xg and canot be reduced to the 
knowledge of the configuration averaged p. The problem is similar to that of the random 
initial condition, and a generalized procedure for (86) depending on the structure of X 
is required. This problem will not be treated presently, as we will limit ourselves to the 
non-random observables. 



1. Averaged one-electron density matrix 



The configuration dependent equation of motion for g can be averaged with the result 
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ih 



[H uv (t),(Q(t))\ = ^ ( ^ R G< 7 G<S A + S<G A - G R S< \ t , 



ft ^(Q(t)) 



ft §^(Q(t)) 



FORW 



ft J^<£?(*)> 



(87) 



and the two terms 



COLL 



The r.h.s. of the Eq. flS8| ) represents the collision rate -^(g(t)) 

correspond to the forward scattering and to the back scattering of the electrons. Altogether, 
a precursor kinetic equation is obtained, from which a true closed kinetic equation for p may 
be derived This equation is not independent of the the equation of motion for G < (and 
its conjugate), but it differs in the direction of integration in the t, t' plane: along the t axis 
for (0), but along the time diagonal t = t' for Eq. flSBp. In a consistent theory, the Green 
functions and the self-energies computed directly should turn (|88|) into a tautology. At this 
point, we may only state that our numerically obtained Green functions obey the precursor 
equation ([58]), postponing all details to a publication about testing the ansatzes. 

From (|8~7D , we see explicitly that in the exact theory the particle number is conserved, 
that is Tr (g)(t) is time independent. A proof will now be given that the SCBA does the 
same, as expected from a conserving approximation, cf. Subsection [IV A 3j . By the Eq. flB8"|) , 
the conserving property is equivalent with the general criterion 



Tr { U R G< - G<E A + E<G A - G R E< } _ # = (89) 

An explicit derivation of the conserving property both for the SCBA and for the CPA was 
previously obtained || in the special case of a rectangular light pulse. Here, we verify that 
SCBA satisfies the Eq. fl39|) for an arbitrary pulse. The SCBA self-energy matrix is given by 
(H). With the propagators written as G R = (G> -G<)$(t-t') and G A = (G < -G > )$(t l -t), 
the expression (|89"D becomes 

t 

( Jdi Tr { X>(t)((G > - G < )V{t)G < - G < V{i){G < - G>) (90) 

V3 < V(t)(G < - G>) - (G > - G < )D(t)G < ) } ) (91) 

The order of the linear operations (•••), / and Tr has been interchanged and the cyclic 
property of trace employed. The Green functions to the left of T>(t) have arguments (t,t), 
those to the right have (i, t) . Clearly, the first term cancels with the fourth term, and the 
second with the third. The result is zero, as required. 



2. Observables and average values 

It will be convenient to compute the averages (|86|) in the Schrodinger picture; this 
amounts to the inverse transformation of the density matrix, 
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^ 0, '° - lnl „ < 92 > 

Qcc ~T Qcv^ i Qvc& ~T Qvv 

In the second line, the expression (0) for O is used. 

The total particle number is a special case of the average value of an observable X = 1 0P , 
which is non-random, so that the averaged p = (g) was sufficient to use in the preceding 
subsection. There is a number of other important observables which are non-random in our 
model. In particular, the electric polarization vector equals to 

P = T 1 iV-iTr(er-<f ?s )) (93) 

because the position vector is represented by the non-random off-diagonal operators r — > 
r cv + r«c- T denotes the volume of the primitive cell. The electric current density 
J = Tq 1 iV _1 er = Tq 1 i\r -1 (i7i) -1 [e:r, H], although a random operator, is also reduced 
to computing (ps>, as 



((j)) = To^-w^-v^s)) 

= T " 1 iV- 1 Tr(er((m)- 1 [H, fe ])) 



(94) 



Ps 



J = l P 



The energy transfer between the external disturbance and the system can similarly be 
simplified in our model in which the optical field has no random component, that is the 
local field corrections are negligible. In that case, the power absorbed per second in a unit 
volume is 

w = E • J (95) 
= To 1 A r_1 Tr(E~^er ((ihy^H, g s })) 

= To 1 N-'Trmr'iH - U(t),H) &)) (96) 

V v ' 

= ^DARK 

We introduce the Hamiltonian 7^ DARK of electrons in the dark sample. The final form of the 
absorbed energy rate is then 

w = r^N- 1 ^- t ((H DARK )) (97) 

The work done on the sample by the light turns into the change of the internal energy of the 
electrons, so that no true dissipation occurs, whether in the exact theory, or in the SCBA. 
This is to be expected in a theory without sources of dissipation. Naturally, this sets an 
upper limit on the time interval in which such purely elastic theory may be valid. 
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B. Computed one-electron properties 



The computed double time Green functions were presented in Sec. [VQ . Here we continue 
with the results of numerical work considering the time-diagonal —ihG < (t, t), that is the one- 
electron density matrix p. 

1. Time evolution of p(t) 

First, we present in Fig. the full density matrix ( |S5|) decomposed into its band- and 
k- vector dependent components. Because in our parabolic band regime the distribution 
is isotropic in the k space, we use the bare electron energy instead of k. In the figure, 
A = e c (k) — e c (k x ), where k x corresponds to the weak field one photon resonance (see 
Fig. 2 and Eq. (|5|)). A is related to the usual detuning by a factor m v /{m c + m v ). For 
the case presented, the spherical layer in the Brillouin zone corresponding to Fig. |T^ is 
given by the inequalities 0.084 < |k| < 0.103 cm., that is 0.134 < \k\/k x < 0.170. Here, 
kx = 2ii/a Q = 0.605 a.u.. 

The diagonal panels show the excitation into the conduction band and the correspond- 
ing depletion of the valence band. The band diagonal elements are real and equal to the 
occupation numbers n CjV (k). The overall effect corresponds to the expectation. A coherent 
transient excites a broad region in the BZ as an energy uncertainty effect. It is followed 
by the persistent excited population in a narrow strip around the resonance, which is the 
strongest for small detunings. 

We use the off-diagonal panels to plot Re p cv and Im p cv which together give both off- 
diagonal elements p cv = p* vc . This is the polarization part of the density matrix. It shows 
a complex oscillatory behavior with an approximate mirror symmetry around the zero de- 
tuning. The polarization is attenuated with the characteristic time h/ImE^ AUK (k), that 
is twice slower than the diagonal transient. This can be understood from fl51|) and fl8"5]) 
by a reasoning similar to that in Sec. |VC2| : the relevant quantities to be compared are 
G^ u (t,t )g(t )G^ c (t ,t) and, say, Gf v (t, t ) g(t )G^ v (t , t) . Both G factors decay in the former 
expression, only one in the latter one. We note that we intentionally do not introduce a 
semi-empirical dephasing time leaving the loss of coherence entirely to the effects of disorder. 

2. Time evolution of observable properties 

Figure [IB] presents, in a matrix arrangement, the total effect of the pulse integrated over 
the BZ, 

pT(t)=N- 1 Y,Pa b (Kt), a,b = c,v (98) 

k 

Full lines in all figures give the result of a direct GF computation. The diagonal panels show 
the total band populations per cell. An almost precise conservation law p™ T + p™ T = 1 is a 
satisfactory result, because the two contributions are computed independently on the basis 
of constituents which do not show much resemblance individually. The coherent parts do 
not compensate (thin line), and the deficit is compensated for by the backscattered flow. 
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Further, we return for the last time to the importance of the induced part of the particle 
self-energy, By a dotted line, we plot the result, if this self-energy component is 

neglected and thus the dark quasiparticle approximation is used. Very clearly, the particle 
number conservation is strongly violated, and the error occurs in the valence band. This is 
in good correspondence with the tests of the semi-group property in Sec. |V B 3| , where the 



most pronounced deviations also occurred for the valence band Green function. This is most 
likely due to the fact that our valence band is not affected by the disorder directly, and its 
behavior is dictated by the indirect influence of the conduction band in a rather sensitive 
manner. 

The off-diagonal component p™ T is shown in the vc corner of Fig. ^ by components, in 
the cv corner in the semi- logarithmic form: p™ T = Re lip . All quantities shown correspond 
to the Galitskii envelope. They are related to the induced polarization per cell by 

T P = ex c ,{p™ T e~ im + cc.} (99) 
= ex cv R(t) cos(f2t + (p(t)) (100) 

Total polarization oscillates rapidly within the envelope ±R. We compare its amplitude 
with the pulse envelope scaled to the rising part of polarization. The pulse and the response 
initially appear to be proportional having constant their ratio and phase difference. Around 
the pulse peak, the response is sub-linear, while the phase changes sign and rises close to hit. 
Polarization vanishes right after the pulse: practically all of the polarization is associated 
with the coherent part of the excitation. The rapid attenuation of the polarization is no 
doubt caused by a cancellation due to dephasing of the rapid oscillations of p vc (k,t) seen in 
the previous figure. 

Using Eqs. (0), (|9"4[ ) and (|100|) , we get for the power absorbed per primitive cell and 



averaged over the Q cycle 

T w = 2Q$(t) {Re(|p uc ) - Q \m Pvc ) 

= 2Q$(t) cos ip) - (2Rsm y?} 

The Q oscillations are rapid compared to variation of the Galitskii amplitudes (RWA condi- 
tion) so that the second term in the braces is dominant. The phase angle between the current 
and the field is then basically if + ~7r. It evolves from a small positive value m 0.0957T to a 
saturation value close to n: the energy is absorbed at first, then partly returned to the field, 
as shown in Fig |T7]. The integral W = /*dt' w(t') multiplied by T is the energy transferred 
between the pulse and one primitive cell up to the time t. It is also shown in the figure. The 
ratio T W (t) / p™ T (t) for times after the pulse gives the energy absorbed per a photoexcited 
electron; it corresponds closely to the basic frequency of the pulse, Q = 1.75 eV. 

Considering the real magnitude of the photoexcitation, we have two electrons per 
primitive cell in the fully occupied valence band, that is the bulk electron density is 
8/al = 3.3 x 10 27 m" 3 . With the fractional excitation ~ 10" 4 , this yields 10 23 10 24 m" 3 
of the electrons excited to the conduction band. Such densities are low enough to justify 
that we neglect the e — e collisions during the time interval considered, cf. @ . The induced 
polarization (|93| ) can be estimated by P = 8/a 3 , P™ T ' ea o as compared with the vacuum 
electric induction D VA c — £oE M - Both quantities are comparable and have the ~ 10 -4 Cm -2 
order of magnitude. 
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VII. CONCLUSION 



In this paper, we have addressed the problem of the response of electrons in semiconduc- 
tors to a short light pulse. Of the three principal scattering mechanisms, we concentrated on 
the impurity scattering, and to make it pronounced, the case of concentrated semiconductor 
alloys was studied. For extremely short times, the disorder scattering in alloys may typically 
be dominant. The scattering regime is easily adjusted by varying the alloy composition. 

The physical picture obtained has the following features: 
oThe alloy scattering in chemically saturated alloys is short range affecting thus large parts 
of the Brillouin zone. 

olt is elastic, which makes some coherence effects pronounced. 

oFinally, it also acts in the dark, so that the dark polaron effect, i.e. use of dressed terminal 
states is essential. Typically, these dark states have a well-defined quasi-particle nature. 
oThe illuminated state cannot be reduced to a redistribution of electrons among the dark 
states, and the induced part of the electron self-energy describes a coherent time non-local 
coupling between the alloy scattering and the excitation process. 

oAs shown by a cumulant analysis, for very short times the dark component describing fast 
quasi-particle formation and the delayed induced component contribute to the self-energy 
additively, so that the formation process is insensitive to the illumination. 
oThe coherence between the light and scattering significantly affects the photoexcited dis- 
tribution; in particular, it conditions the particle number and energy conservation. 

All these conclusions are obtained analytically and confirmed numerically. 

The paper discusses also technical questions. The NGF describing the process are ob- 
tained by a direct solution of equations equivalent to the Kadanoff-Baym equations. The 
asymmetric LW choice of computing first the propagators and using these as an input in the 
equation for the particle correlation function is advantageous in the present case of an elastic 
scattering. The isotropic scattering would require /c-integration over the whole BZ. This is 
overcome by separating out the dark self-energy and integrating explicitly only the induced 
part, stemming from a narrow slice around the resonance in the /c-space. To suppress rapid 
oscillations of the integrand, we work in the Galitskii representation, in which the pulse is 
represented by its envelope, while the high basic frequency enters only as a relative shift of 
the bands. 

For the future work, this paper suggests several directions as an immediate sequel: 

(i) It will be important to extend the NGF approach to the Coherent potential approxima- 
tion, which is a self-consistent intermediate theory suitable for all realistic alloy scattering 
strengths. 

(ii) Then random ("correlated") initial conditions will play a non-trivial role, and should be 
incorporated. 

(iii) These analytical approximations for the NGF may be, in the alloy case, compared with 
a direct simulation for explicitly generated random alloy configurations and a numerically 
performed configuration average. 

Independently of these extensions and modifications, the results of the present work 
can readily be used for an analysis of various methods to develop quantum kinetic equations 
based on an ansatz generalizing the KBA, and for a comparison of the resulting approximate 
solutions with the NGF results. Such analysis should serve to find the validity limits of the 
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individual ansatzes and the underlying physical explanation. A communication concerning 
these questions is also in preparation. 
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APPENDIX A: G< WITH INITIAL CORRELATIONS 



Like in the main text, the basic definition of G < is the Eq. (52), but now with the 
initial condition (|53|)at the initial moment to being random. We express the propagators as 
QR = G R + G R {V-E R )g R , g A = G A + G A (V - E A )G A . The expression G< = ih{g R oPo g A ) 
becomes 



G< = ih [G R ( p )G A 

+ G^U^X? - E A ))G A + G R ((V - s R )g R tPo )G A 
+ G R ((V - E R )g R op o g A {V - E A ))G A ) (Al) 
it > t , t' > t ) (A2) 



The three lines of ( |A^ ) correspond to an averaged (" uncorrelated" ) initial condition, to 
the random field analogue of the correlated initial conditions, and to the genuine particle 
correlation, respectively. These terms have the same general structure G R . . . G A , but they 
differ in the number of inner times equal to to, two, one, and none. We make this explicit 
using the " " subscripts and define 

D A< = iHop^iV - E A )) (A3) 

A< ee ih((V - Z^cpo) (A4) 
E< ee {{V - E R )g<{V - E A )) (A5) 

With these definitions, the final form of G K is 

G< = ihG R { c p )G A + G R C A < G A + G R k<G A + G R E < G A (A6) 



This equation extends the Eq. (p4|) by incorporating a memory effect caused by the disorder 
induced initial correlations. In all other respects, comments made to (p4"D apply here equally 
well. 

The diferential equation analogous to the Dyson equation (|6"4]) is obtained, if we differ- 
entiate {"KB) from the left, and employ (^l|) for G R in the form 



(ihd t -H MF )G R = E R G R + I. 
For t > to, the 5-singularity is effective for only two terms of ( |A6|) , and the resulting equation 



reads 
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(ihd t - iJ MF )G < = E R G < + E<G A + A<G 



(A7) 



The two seemingly missing terms reappear in the initial condition: 



G<(t ,t') = (oPo)G A + ( A < G A 



t=to,t' 



(A8) 



It remains to develop an approximation scheme for the self-energy, both for the regular part 
E < and for the singular part D A < + A< . The two A's are equivalent: 



A < (t,t , ) = -{A<(t',t)y 



(A9) 



so that the whole < component of the self-energy obeys the symmetry (fig). Considering, 
say, Aq , we see from (K4) that it only involves retarded quantities, which do not depend 
on the initial conditions for elastic scattering. Thus, A< is a linear functional of the initial 
distribution p , and it does not involve any particle- hole correlation. In other words, A's 
can computed on the propagator level. The averaging procedure leading to ([A4D depends 
on the initial distribution, however. 



The Eq. (A5) for S < does not explicitly depend on the initial condition. In the SCBA, it 
will have the form £< BA = {VG<V) identical with (||), but with G< given by ([£§) and (|AJ. 
Only in higher order approximations, the genuine correlations of the three-terminal type, 
like (T>G R (op — (op ))G A V) , start playing their role, and must be systematically accounted 
for. 



APPENDIX B: RANDOM OBSERVABLES 

Let us return to but write it in the Schrodinger picture, cf. (P2|): 

((X)) = Tr(Xg s } (@3) 

If the observable X is random itself, the double average has to be performed at the end, as 
discussed in the main text. The average can be given two different forms: 

((X)) t = h 2 Tr(xg R (t, i )^o£ A (to, *)> (Bl) 
= -ihTr(Xg < (t,t' -> t)) (B2) 

The first form, Eq. ([Bl]), makes explicit that the random initial condition and a random 
observable enter the averages in a symmetrical way, and any general technique for their 
evaluation has to be built up symmetrically. The second form of ((X)), seemingly trivial, 
replaces the density matrix by the particle correlation function, and this may lead to explicit 
averaging procedures in special cases. 

As an important example, let us consider the dark Hamiltonian in our model and calcu- 
late pi dark))) a quantity related to the Joule work by (|97|). This does not go with (p), but 
we may proceed as follows: 
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DARK 

+ u(t))g< 

m|G< - U{t)G< 

Tr(n BARK (-in)g < (t,t' 

ft 2 Tr(|G<(M'^£)) 



('^DARK^ < ) 

((T^dark))* 



Tr(f/(i)p) 



(B3) 



Use is made of the fact that U is non-random, and its average reduces to the use of p = {g). 
The essential point is, however, that the derivative on the r.h.s. is NOT a derivative of 
the density matrix. The trick is similar to the computation of correlation energy using the 
one-particle GF in usual many-body theory The derivative can be found directly (it would 
be enough to store it during the solution of the Dyson eq. for G < ), or it could be expressed 
using this DE (see Eq. 

Thus, we can do two things: either get ((H DARK ))t directly from (|E3|), or express it in terms 
of the integral of the " Joule heat" , W{t) = fdt' w{?) , using (|97|) . The two results should be 
the same. This provides another criterion for the conserving nature of any approximations 
involved in the calculation. 
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TABLES 



TABLE I. Concentration dependence of the alloy electron structure: weak field one photon 
resonant energy , level broadening and particle lifetime 



c 




hCl x /eV 


|Imz fc |/ eV 


r/fs 


0.00 




1.80 


0.0000 


oo 


0.05 




1.75 


0.0020 


165 


0.15 




1.63 


0.0059 


53 


0.50 




1.34 


0.0105 


31 


TABLE II. Lowest moments and cumulants: exact and in the SCBA 


V 


M p SCBA 


M p EXACT 


C p SCBA Cp EXACT 





1 


1 


1 1 




1 













2 


7 


7 


7 7 




3 


7A»i 


7/ii + 7(c B - c A )6 


7Mi 7/"i + 7(cb ■ 


- c A )5 


4 


7/x 2 + 27 2 


7^2 + 27(c B - c A )5fii+ 


7^ 2 -l7 2 7A*2 + 27(c B - 


c A )5[xi+ 






+ 7 (5 2 -3 7 ) 




-67) 


TABLE III. Characteristic parameters of the pulse 


t P [ps] 




Q[eV] E M [Vm _1 ] 


^FWHM [PS] 




0.1 




0.01 1.89xl0 7 


0.121 


1.93^ 
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FIGURES 



FIG. 1. Local Green function F Q (E + iO) of the pure A crystal for energies between the band 
edges 1.5 eV and 13.5 eV on the real axis. 

FIG. 2. Renormalized dispersion law function of bare energy e/t for four values of 

concentration c = 0.; 0.05; 0.17; 0.5. Dotted lines: mean field rigid shift of the whole band. Thick 
line: full dispersion law with polaron shift included. Energy dependent band broadening shown by 
lining each of the quasi-particle energy by thin line at a distance zblmz^. Dot-dashed lines: shifted 
valence band e„(k) + ?i$7 x (c). 

FIG. 3. Renormalization constant as a function of bare energy for the same values of concen- 
tration as in Fig. 2. Upper part (thick lines): modulus \Zk\ — 1. Lower part (thin lines): phase 
arg(2' fc ). 

FIG. 4. Evolution of the "action" s(t) for dark propagator on femtosecond time scale. Fig. 4a: 
real part, Fig. 4b: imaginary part. Full line: full SCBA solution. Dashed line: pole approximation 
( |30| ) in time domain. Detail for the shortest times < 0.3 fs (small windows) expanded in the insert 
in the left lower corner and compared with the lowest (full dots) and next lowest (open dots) order 
cumulant expansions. The insert in the right upper corner: effective self-energy <jfc, Eq. (^) (full 
line) tends rapidly to the constant quasi-particle (pole) energy (dashed line). 

FIG. 5. Schematic flow chart for calculation of T,f(t,t'). Input: dark QP for p = 0. Iteration 
counter: p. Output: Tif(t,t') and G R (k;t, t') (not indicated). 

FIG. 6. Schematic flow-chart for calculation of U < (t,t'). Input: propagator components U R , 
^ind> G A . Output: E < (t,t') (and G < (k;t, t'), if required - not shown). Main differences with 
respect to E R D (t, t') (see Fig. 5): o non-iterative process (see text); o the initial time t' . . . outermost 
loop, the current time t . . . innermost loop. This "normal" order of loops was reverted for S R . 

FIG. 7. Imaginary part of the induced self-energy (^) as a function of both time variables. 
Insert: the same plot rotated to the view in the direction of the arrow (along the t + t' axis). 
Vertical scale unit: eV/ps. 

FIG. 8. Band diagonal elements |i?iG^| 2 for resonant wave vector k x as a function of t and t'. 
Upper panel: b = c. Lower panel: b = v. In each panel, two surfaces are marked by dots, the 
full SCBA G R , and its approximation with E R neglected. As an aid for eye, the exact surface has 
several lines t' = const marked by heavier points. The dot patterns permit to view the surfaces 
either in the t, t' coordinates, or in the Wigner coordinates t ± t' . 
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FIG. 9. Test of the semigroup multiplicative property (|83|), Heavy line: iHG R (t, — 0.2). Thin 
lines: \TiG R (t, t")ihG R (t", -0.2) for t > t". Quantity plotted: \ihG cc \ 2 + \ihG vc \ 2 (upper panels) 
and \ihG vv \ 2 + \ihG cv \ 2 (lower panels). The pulse envelope (||) and the factorization times t" 
(squares) shown at the top. The first, third, and fifth t" were singled out also in Fig. ||. Left 
hand column: full SCBA, Sj n ^(t,t') included. Right hand column: WW approximation, Ej n( j(i, i') 
neglected. Here, the factorization is valid and the thin lines merge with the basic heavy plot. 



FIG. 10. Self-energy £ < (t,t') as a self-consistent result of eqs. (|72j) and (|77|). Upper panel: 
real part, lower panel: imaginary part. Inserts: the same plot viewed in the direction of the arrow 
(along the t + t' axis). Units of E: eV/ps. 

FIG. 11. The particle correlation function G^.(k;t,t') for resonant detuning k = k x . Upper 
panel: real part, lower panel: imaginary part. Inserts: view along the arrow - as in FigJlO|. Units: 
fiG is dimensionless 

FIG. 12. lm.Gf c (k x ;t,t') of the previous figure split into its coherent (upper panel) and inco- 
herent (lower panel) parts according to (|54|). 

FIG. 13. The coherent part ImG< C0H (k x ; t, t') of Fig. [L2| (upper surface, coarse grid) contrasted 
with the calculation, where U R D (t,t') was neglected. (lower surface, fine grid). 



FIG. 14. Time diagonal ImG^,(k x ; t, t' = t). Thick line: data of Fig. 11. Thin line: the same 
quantity, but 2J R D neglected in computation. 

FIG. 15. One-electron density matrix decomposed into its band and A:-vector dependent com- 
ponents. Detuning A equals to the excess energy e c (k) — e c (k x ) measured with respect to the 
one photon resonance. Normalization to one primitive cell, p is dimensionless. The A = profile 



identical with the thick line of Fig. 14 



FIG. 16. Total photoexcitation p^ T , Eq. (^). Diagonal panels: total band population per cell 
(thick line ... direct GF computation; thin line ... coherent part only; dashed line ... computation 
neglecting E R D ). Conservation law / o™ T + p™ T = 1: very good for direct GF computation, poor 
in both other cases. Off-diagonal panels: The vc corner: the off-diagonal component p™ T . The cv 
corner: the same quantity in semi-logarithmic form /9™ T = -Re 1</3 . Thick line . . . R. Thin line . . . (p. 
Dotted line . . . pulse envelope $ scaled to coincide with R at early times. Note the good fit. All 
quantities shown are dimensionless. 

FIG. 17. Power absorbed per primitive cell and averaged over the f2 cycle ... thick line; integral 
energy transfer between the pulse and one primitive cell ... thin line Units: eV/ps for T w, eV for 
T W. 
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